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Scanning tunneling spectroscopy (STS) is a useful probe for studying the cuprates in the super- 
conducting and pseudogap states. Here we present a theoretical study of the Z-map, defined as 
the ratio of the local density of states at positive and negative bias energies, which frequently is 
used to analyze STS data. We show how the evolution of the quasiparticle interference peaks in the 
Fourier transform Z-map can be understood by considering different types of impurity scatterers, 
as well as particle-hole asymmetry in the underlying bandstructure. We also explore the effects of 
density wave orders, and show that the Fourier transform Z-map may be used to both detect and 
distinguish between them. 

PACS numbers: 74.55,+v, 74.72.-h, 71.45.Lr 



I. INTRODUCTION 

Conventional metals and superconductors exhibit 
particle- hole symmetry at low energies, which is well de- 
scribed by Fermi liquid and BCS theory. In contrast, a 
notable feature of the high-T c cuprate superconductors 
is the number of ways in which they break particle-hole 
symmetry. Electron or hole doping the Mott insulating 
parent compounds leads to superconductivity, although 
the maximum T c is drastically different, reaching 150 K 
for hole doped compounds, while only 30 K in the case of 
electron doping. Another example is a pronounced asym- 
metry in the conductance spectrum for tunneling elec- 
trons and holes, measured by scanning tunneling spec- 
troscopy (STS). 1 These asymmetries can be understood 
in the context of doping a Mott insulator, where the ad- 
dition of electrons is more difficult than extraction, due 
to the strong electron-electron repulsionP^ 

Recent experiments have observed further particle- 
hole symmetry breaking in the pseudogap state, a phase 
whose origin is debated, but may be due to preformed 
Cooper pairs lacking phase coherence, or a compet- 
ing order. 4 Angle-resolved photoemission spectroscopy 
(ARPES) 5 in the pseudogap regime reveals particle-hole 
symmetry breaking in the form of a bend-back in the 
energy dispersion misaligned with the Fermi momen- 
tum ki?, which can arise from a density wave (DW) 
order. 6 7 Because superconducting fluctuations should 
be roughly particle-hole symmetric, while spin and 
charge orders generally induce particle-hole asymmetric 
changes, searching for particle- hole symmetry breaking 
is a useful way to distinguish between superconductivity 
and DW orders. 

In this paper we explore the signatures of particle-hole 
asymmetry in the local density of states (LDOS) mea- 
sured by STS, a technique which has been used exten- 
sively to study the cuprates in both the superconduct- 



ing and pseudogap phases PdaCD STS has revealed that 
the cuprates have a spatially inhomogeneous electronic 
structure, including modulations in the LDOS and su- 
perconducting gap magnitude GSHUl i n the d-wave super- 
conducting phase, the LDOS modulations can arise from 
quasiparticle interference (QPI), due to the s cattering 
of wave-like quasiparticles off impurities! 12 * 15 * 22 -^^ The 
wavevectors of the modulations can be determined from 
the Fourier transform of the LDOS. 

Fourier transform STS (FT-STS) data are often ana- 
lyzed using the Z-map, defined as the ratio of the LDOS 
at positive and negative bias energies, 



n(r, +w) 
n(r, — lj) ' 



(i) 



Experimentally, the Z-map has been used to cancel both 
non-dispersing "checkerboard" modulations and system- 
atic errors due to tip elevation uncert ainty, in order to re- 
veal the superconducting QPI pattern!^ 2 ^ By definition, 
the magnitude of the Z-map clearly carries information 
on particle- hole asymmetry in the LDOS. 

Due to superconducting coherence factors, the Z-map 
in fact enhances the intensity of LDOS modulations from 
QPI, as discussed in Ref. |28j To summarize their argu- 
ment, since STS tunnels electrons rather than quasipar- 
ticles, the strength of the measured conductance at r 
depends on the magnitude of the hole and electron am- 
plitudes, |ix n (r)| 2 and |v n (r)| 2 (n labels the eigenvalue of 
excitations). Expressing the LDOS in terms of the elec- 
tron Green's function n(r,cj) = — (l/7r)JraG?ii(r, lj), the 
dependence of the LDOS on coherence factors is clear: 

n(r,w>0)=J2\u n (r)\ 2 S(w-E n ) (2) 

n 

n(r,Lj < 0) = \v n (r)\ 2 S(LU + 
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Assuming well denned excitations in energy (so the sum 
over n includes just one term), the Z-map reduces to 



By recalling the sum rule J2 n \ u n( r )\ 2 + l^n(r)| 2 = 1, it 
is clear that whenever |^(r)| 2 is large |^(r)| 2 is small, so 
taking the ratio of the coherence factors will enhance the 
maxima. 

The superconducting QPI pattern in the cuprates is 
dominated b y scat tering wavevectors predicted by the 
octet modelP3S3 The octet model is derived by not- 
ing that for a d-wave superconductor, the quasiparti- 
cle dispersion exhibits banana-shaped contours of con- 
stant energy (CCE). The momentum space LDOS is 
maximal at the tips of the CCE, where the curvature 
is strongest. Quasiparticle scattering is dominated by 
wavevectors qi — q^j that connect the endpoints of the 
CCE, depicted in Fig. [I] The presence of particle-hole 
symmetric wavevectors dispersing according to the octet 
model is evidence that an observed gap is due to super- 
conductivity; the shape of the underlying Fermi surface 
and the momentum dependence of the gap can be deter- 
mined from this dispersion. 

The evolution of the QPI intensity with bias energy 
has recently been a subject of discussion in the litera- 
ture. FT-STS experiments have observed that the in- 
tensity at some octet q- vectors disappears when the tips 
of the banana-shaped CCE cross the antiferromagnetic 
zone boundary (AFZB), which has been interpreted as a 
signature of the loss of QPI, or more generally, the loss 
of quasiparticles) 2 ^ In contrast, recent ARPES measure- 
ments observed well defined quasiparticles all the way 
from the node to the antinode, and it has been argued 
that the extinction of QPI at the AFZB observed in FT- 
STS can be attributed simply to the momentum depen- 
dence of impurity scattering P"J Additionally, the presence 
of a competing spin DW order coexisting with supercon- 
ductivity can explain the disappearance of QPI beyond 
theAFZBP 

In addition to dispersing octet q-vectors from su- 
perconducting QPI, static quasi-periodic "checkerboard" 
modulations have been observe d in both the supe rcon- 
ducting and pseudogap phases. ™E2HSl2l^^ A1 _ 

though the origin of these modulations remains unknown, 
there has been much discussion about their relation to 
a possible DW order in the pseudogap phase; scenarios 
that have been proposed include orbital current induced 
d-density waves pa one dimensional stripes,^ nematic or- 
der,^ short range charge order connected to nested p arts 
of the Fermi surface in antinodal regions IMlMEil anc [ 
disorder-induced charge orders p2 

We present a detailed study of what can be learned 
from the Z-map in superconducting, DW, and coexisting 
phases. In the first part of the paper we show that the 
evolution of the Z-map intensity with bias energy and 
doping reflects the underlying particle-hole asymmetry 




FIG. 1. The octet model. Scattering is dominated by 
wavevectors qi-q7 which connect the eight points at the ends 
of banana-shaped contours of constat energy (red curves) 
centered on the normal state Fermi surface (dashed line). 
The phase of the d-wave superconducting order parameter 
is shown in green (positive) and orange (negative). 



of the bandstructure. We also compare the QPI pat- 
terns from different impurity types, and show how the 
Z-map can distinguish between them. In the second half 
of the paper, we consider DW order- both long range and 
fluctuating- and show that due to inherent particle-hole 
symmetry breaking, DW order produces a unique pattern 
in the Z-map. As a result, the Z-map may be a useful in 
detecting and distinguishing between DW orders. 

The plan for the rest of this paper is the following: 
in Section II we discuss general properties of the Z-map 
and its connection to particle-hole asymmetry. In Section 
III, using a single site impurity, we study the evolution 
of the QPI intensity with bias energy and doping, while 
in Section IV, we introduce an extended impurity with 
scattering treated self consistently. We investigate the 
signatures of DW order in the Z-map in Sections V-VII, 
considering the effect of long range and fluctuating DW 
orders, and their coexistence with superconductivity. We 
conclude in Section VIII. 



II. PARTICLE-HOLE ASYMMETRY IN THE 
Z-MAP 

Because the CCE in the superconducting state are 
particle-hole symmetric, the octet wavevectors are the 
same at positive and negative bias energy (q^(cj) = 
q_i(— a;)), however, the QPI intensity (the magnitude of 
the LDOS modulations) at those wavevectors need not 
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be. Here, we consider the effect of inherent particle-hole 
asymmetry in the underlying cuprate bandstructure on 
the Z-map intensity; the difference between injecting and 
extracting electrons from a Mott insulator would further 
enhance this asymmetry in a real material. The Z-map 
intensity evolves with bias energy in part because n(r, lj) 
is closer to particle-hole symmetric at low energies than 
at higher energies (where there is more asymmetry due 
to the proximity of the van Hove singularity). Addition- 
ally, variation in the QPI intensity with doping can be 
understood by appealing to particle-hole asymmetry, as 
doping moves the van Hove singularity. 

In order to quantify these ideas, we first decompose 
the Z-map into even and odd functions of bias energy, 
which makes clear the dependence on the particle-hole 
asymmetry of the underlying bandstructure. The LDOS 
can be written as the sum of a uniform part and a part 
modulated by impurity scattering: n(r,co>) = no(uj) + 
Sn(r,(j). When LDOS modulations are weak, such that 
\Sn(r : Lj)\ « | no (a;) |, the Fourier transform Z-map can 
be expanded as 4 ^ 



£n(q, lj) (5n(q, — cj) n 
n (uj) n Q (-u) * 



= Z (lj)(S^ + 5Z(q,(j)) 

where Zq(uj) = tiq{lj) / uq{—lj) . Decomposing the 
LDOS as £n(q, lj) = 5n even (q,oj) + 5n odd (q, where 
Sn even (q,uj) = Sn even (q, -uj) and Sn odd (q, lj) = 
—5nodd(q., — w), the Z-map modulation can be expressed 
as 

SZ(q,(j) = a(u)5n odd (<i,uj) + j3(uj)5n even (^uj) (4) 



where 



and 



a(uj) 



1 



i 



n (cj) n (-Lj) 



n (uj) no(-u)' 



(5) 



(6) 



The factors a(uj) and /3(lj) reflect the amount of 
particle-hole asymmetry in the underlying band struc- 
ture. For a particle-hole symmetric bandstructure, 
a(uj) = 2/no(uj) and /3(lj) = 0, so the Z-map modula- 
tion depends only on 5n odd (q, lj): 



n (uj) 



(7) 



The cuprate bandstructure weakly breaks particle-hole 
symmetry for energies less than the gap maximum (lj < 
A ), so in this case a(u) » (3(lj). As a result, the 
Z-map modulation is dominated by Sn odd (^Lj) in the 
cuprates. We use the expressions derived in this section 
frequently in our subsequent discussion. 



III. SINGLE SITE IMPURITY SCATTERING 

Since QPI requires the presence of impurity scatter- 
ers, in this section we lay out the formalism for single 
impurity scattering, while in later sections we treat an 
extended impurity model. We consider three impurity 
types in this paper: those that modulate the d-wave 
superconducting gap, the bond (nearest neighbor hop- 
ping) parameter, and the site energy. The QPI pattern 
observed in FT-STS experiments arises from scattering 
from multiple impurity types, which can be difficult to 
disentangle. However, recent experimental and theo- 
retical 40 work studying the effect of a magnetic field on 
QPI has shown that vortices induce gap inhomogeneities 
which can scatter quasiparticles. The magnetic field does 
not create additional scattering vectors, but enhances in- 
tensity at octet q-vectors which originates from scatter- 
ing from gap inhomogeneities, so these wavevectors now 
can be isolated and studied in detail. 

The theory of QPI from various impurity types has 
been studied by several authors. Cal culation s of QPI 
from single site and multiple impurities l 29 l 32 l 42 t El repro- 
duce the locations of experimentally measured wavevec- 
tors, but the features in the QPI pattern are "streak-like" 
rather than "point-like" as in experiment. Zhu et.al ^ 
studied QPI from multiple impurities, and highlighted 
the need for more realistic impurity models. Other works 
have shown that many of the discrepancies between the- 
ory and experiment can be resolved by considering an 
extended long range impurity model, 46 or by including 
a spatially inhomogeneous chemical potential and super- 
conducting gapP^ However, the aim of this paper is to 
focus on the generic behavior of the Z-map, and we ne- 
glect more realistic impurity modeling. 

We now review the formalism for single site impurity 
QPI calculations. While we neglect electronic correla- 
tions in this work, we do not expect them to change our 
overall conclusions, which largely are derived in a fairly 
straightforward way based on symmetry arguments. In 
addition, we confine our considerations to the low tem- 
perature phases where there is strong experimental ev- 
idence that there are well-defined quasiparticles. The 
BCS Hamiltonian in the Nambu formalism is given by 



(8) 



where ip^ = (cL.,c_k^) and 



A k 



A k -e k 



Here c£ creates an electron with momentum k and spin 
a, ek is the bandstructure, and for our considerations 
of the cuprates, we use the d-wave superconducting gap 
Ak = Ao(cos/c x — cosk y )/2. The Green's function is 
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given by 

G(k,iw„) = (ituj - Ak)" 1 (9) 
_ ioJnTp + A k fi + e k T 3 
(iu n )i - El ' 

where fo is the identity matrix, fi } 3 are Pauli matrices, 
and£ 2 =e 2 +A 2 . 

On a lattice of N sites, scattering from a single point- 
like impurity at site (0,0) is described by the Hamiltonian 



Himp = ^2 ^^,0^0 + h.C. 



(10) 



where ipl = (cj^, c r ±) and T r = T\^T r Q is the scattering 
T-matrix. Contributions to the T-matrix from different 
impurity types are classified by how they modify elec- 
tron parameters: bond modulations St and site-energy 
modulations S/i occur in the T3 channel, while d-wave 
superconducting gap modulations 5 A occurs in the fi 
channel. 

Superconducting gap and bond modulations modify 
the parameters on the four bonds immediately surround- 
ing the impurity site. On a lattice with spacing a the 
real space T-matrix is given by 

T r5 o = 5m[5(r-ax)±5(r-ay)+5(r+ax)±5(r+ay)} (11) 

where 5m = St(SA) is the amplitude of the modulation 
and the +(-) is for bond (<i-wave gap) modulation. Site- 
energy modulation occurs only at the impurity site, for 
which T r o = SfjiS(r — 0). Transforming to momentum 
space, the T-matrix is 



T£, k+q = *Ari(A k + A k+q ) 
r k,k+q = Stf 3 (t k + *k+q) 



[ k,k+q 



S/lf 3 



(12) 

(13) 
(14) 



for gap, bond, and site energy modulations, respectively. 
The definition of JA k was given above, and t k = cos k x + 
cos k y . 

From the momentum space form of the T-matrix, 
we can already make two important observations about 
the energy dependence of the QPI intensity. 31 First, 
T£ k+q = for all octet q- vectors when k and k + q lie 
along the AFZB. Therefore, for impurities that modu- 
late the bond parameter, extinction of QPI at the AFZB 
is simply due to the momentum-space form of the T- 
matrix, rather than a loss of quasiparticles. Including 
next-nearest-neighbor bond modulations lifts the com- 
plete extinction of QPI intensity at the AFZB, although 
these modulations have a smaller amplitude. Addition- 
ally, for superconducting gap modulations, T^ 



k,k+q 







for octet q-vectors that connect points of opposite gap 
phase (q2, q3, q6, qz). All wavevectors are present for 
site energy modulations. 



Born scattering in a particle-hole symmetric 
band 



We now classify the different impurity types by 
whether they create an LDOS modulation that is even 
or odd with respect to bias energy, for the simple case of 
scattering in a particle-hole symmetric band. This simple 
classification scheme will help us understand the behav- 
ior of £Z(q, uj) in the cuprates, where this classification 
still holds, but only approximately. The symmetry of the 
LDOS modulation due to superconducting gap and site 
energy variations already was described in Ref. 48 using 
the Bogoliubov-de Gennes equations for a <i-wave super- 
conductor with classical phase fluctuations. We take a 
different approach: we derive the LDOS modulation due 
to scattering from a single impurity site using a momen- 
tum space T-matrix formalism. 

The LDOS modulation for scattering by wavevector q 



c 30 



£n(q, uj) 



d 2 k 

(2tt)2 

Im[G(k, lo 



- i£)Xk,k+qG(k + q, w 



(15) 
■iS)]u 



where G(k, uj + i5) is the analytic continuation of the 
Nambu Green's function to real frequencies. This ex- 
pression can be separated into parts that are even and 
odd with respect to bias energy uj: 



Jn(q, uj) 



d 2 k S(uj - £ k ) + 5(uj + E k ) 



(2tt)2 

£ra[/i(k,q) + / 2 (k, q)/uj 



(16) 



where £m, /i(k, q) and ^(k, q) depend on the type of 
scatterer. For gap modulation 

/f(k,q) = (A k + A k+q ) 2 (17) 
/f (k, q) = (A k + A k+q )(A k e k+q + A k+q e k ), 

while for bond modulation 

A b (k,q) = (t k + t k+q )(e k + e k+q ) (18) 
f\ (k, q) = (t k + £ k+q )(cj 2 + e k e k+q - A k A k+q ), 

and for site energy modulation 

A s (k,q)=e k + e k+q (19) 
fi (k, q) = uj 2 + e k e k+q - A k A k+q . 

Noting that for a particle-hole symmetric bandstructure, 
e k+Q = — e k for Q = (7r,7r), it can be shown that 



/!(k,q)=0 
/ 2 b (k,q)=0 
/l S (k,q)=0 



(20) 



so that for superconducting gap and bond modulations, 
the LDOS modulation is purely even with respect to bias 
energy, while for site energy modulation, the LDOS mod- 
ulation is purely odd. 
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As a result, since 5Z(c[, uj) depends only on £n oc ^(q, uj) 
for a particle-hole symmetric bandstructure, only site en- 
ergy modulations contribute to the spatially varying Z- 
map. When p-h symmetry is broken, all three types of 
impurities cause both odd and even modulations in the 
LDOS. 



Born scattering in a particle-hole asymmetric 
band 



Having classified the LDOS modulations in a particle- 
hole symmetric band, we now turn to a particle-hole 
asymmetric bandstructure describing the cuprates: 

ek = — 2t(cos k x + cos k y ) — At' cos k x cos k y 
— 2t // (cos 2k x + cos 2k y ) — 

4.111/ ^ I I 7 



— At (COS ZK X + COS ZKy) — 

At'" (cos 2k x cos k y + cos k x cos 2k y ) — ji 



(21) 



where t, t', t" , H" , = 0.22, -0.034315, 0.035977, 
0.0071637 eV respectively and the chemical potential \i is 
adjusted to control the filling. This bandstructure is ob- 
tained from a tight binding fit to ARPES data on nearly 
optimally doped Pb-Bi2201 described in Ref. [71 The d- 
wave superconducting gap maximum is set to A =35 
meV. The uniform DOS no(uj) for this band structure 
in the superconducting state is shown in Figure [2| a). 

Fig. P[b) shows the factors a(uj) and /3(uj) (defined 
in Eqs!~[5] and [6|, which control the contributions of 
5n dd(q,a;) and Sn even (q, uj) to 5Z(q,u;), respectively. 
The term a(uj) is an order of magnitude larger than 
f3(uj), strongly enhancing the contribution from Sn dd- 
Since no (us) increases with energy |cj|, especially at neg- 
ative bias due to proximity to the van Hove singular- 
ity, a(uj) decreases. Thus the Z-map intensity should 
decrease with increasing energy, unless the strength of 
^^odd(q, u) increases markedly. 

Figure |3] shows \5n dd(q., <^)|, as a function of scattering 
wavevector q and energy uj for superconducting gap (pan- 
els a-b), bond (c-d), and site energy (e-f) modulations. 
Two Brillouin zone cuts are shown, which highlight dif- 
ferent octet q- vectors (refer to Fig. [I] for the definitions of 
these vectors): qi and qs in an anti- nodal (0, 0) — » (7r, 0) 
cut (panels a, c, and e), and q3 and q7 along a nodal 
(0,0) — » (7r,7r) cut (panels b, d, and f). 

In the anti- nodal cut, qi and qs start out with negli- 
gible intensity at low energy, and then increase in inten- 
sity as energy increases for all three types of modulation 
(panels a, c, and e). The effect is most pronounced for 
superconducting gap modulation, where the intensity in- 
creases by three orders of magnitude. For the nodal cut, 
wavevectors q3 and q^j are not present for superconduct- 
ing gap modulation (panel b), due to the momentum de- 
pendence discussed at the beginning of this section. For 
bond modulation (panel d), q3 and q7 both start out 
with low intensity at small energies, and then quickly 
grow in intensity. This intensity eventually starts to de- 
crease with increasing energy. The behavior of q3 and q7 
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FIG. 2. (a) Density of states for the band structure defined 
in Eq. |21| (b) ol(uj) and /3(cj) defined in Eqs. [5] and [6] for 
< n >=0.73. 



in the case of site energy modulation (panel f) is qual- 
itatively similar, although both wavevectors continue to 
gain intensity at high energies. 

In summary, from simple DOS arguments, the Z-map 
intensity 5Z(q_, uj) at all octet wavevectors should de- 
crease with increasing |cj| due to the reduction in a(uj). 
However, a strong increase in (5n ^(q, a;), for example, 
at large \uj\ in the case of gap modulation, or a peak at 
intermediate |cj| in the case of bond and site energy mod- 
ulations, could outweigh this behavior, depending on the 
overall strength of the impurity scattering. The overall 
evolution of the Z-map intensity will depend on the in- 
terplay of these two effects, which we will illustrate in the 
next section. 



IV. SCATTERING FROM AN IMPURITY 
PATCH IN SELF-CONSISTENT T-MATRIX 
APPROACH 

To more accurately describe impurity scattering and 
better compare to FT-STS data, we consider the effect 



6 




FIG. 3. Odd part of the LDOS modulation | Sn dd (q, u)\ 
as a function of scattering wavevector q and energy uj for (a, 
b) superconducting gap modulation, (c, d) bond modulation, 
and (e, f) site energy modulation. The wavevector q for the 
left column (panels a, c, and e) is along the anti-nodal cut 
(Q — (0, 0) — >> (7r,0)), while in the right column (panels b, d, 
f) it lies along the nodal cut (q = (0, 0) — > (tt, 7r)). The boxes 
mark the location of the q-vectors as calculated in the octet 
model that lie along the antinodal (qi and qs) or nodal cuts 
(q3 and qr). Intensities are in arbitrary units. 



of an extended impurity! 31 * 49 * We use an extended patch, 
which is embedded in a finite periodic lattice, and cal- 
culate the T-matrix self-consistently in real space. The 
real space Green's function is given by 



G(r,iuj n ) 



(22) 



where G(k, iuj n ) is the Nambu Green's function given in 
Eq.jiJ The initial real space T-matrix Tb(ri, r 2 ) gives the 
amplitude for scattering between two adjacent sites in 
the patch for superconducting gap or bond modulation, 
or the amplitude for on-site scattering for site energy 
modulation. The T-matrix for scattering between any 
points i*i and T2 in the impurity patch is then determined 
self-consistently from 

f(n, r 2 ) = Tb(n, r 2 ) + Tb(ri, r')G(r' - r)f (r, r 2 ) 

(23) 



(the uj dependence of the Green's function has been sup- 
pressed). The LDOS modulation in real space due to 
quasiparticle scattering is given by 

Sn(r,u) = --ImJ2 (g(t - n)f (n, r 2 )G(r 2 - r)) 

7T L — ' V /ll 

ri,r 2 

(24) 

and the real space Z-map is 

no(u) + 5n(r, uj) 



Z(t : uj) = 



no(-uj) + (5n(r, — uj) 
The Fourier transform (FT) Z-map is then 

Z(q, W ) = l^e-^ r Z(r, W ). 



(25) 



(26) 



The figures shown in the rest of this paper were made 
using a 9 x 9 extended impurity embedded in the cen- 
ter of a 512 x 512 real space lattice. The magnitude of 
the modulation parameters in the initial T-matrix are a 
maximum on the bonds in the center of the patch, and 
decrease in magnitude within a "Gaussian-like" envelope, 
in order to minimize the effect of "ringing" that would 
occur from an abrupt fall-off of the T-matrix at the patch 
boundaries. Since the impurity in our calculation is em- 
bedded in the center of the real space lattice, it creates 
an LDOS modulation Jn(r,cj), and hence a real space Z- 
map Z(r, uj) that is C4 symmetric. As can be seen from 
Eq.|26j the Fourier transform of a C4 symmetric quantity 
(here Zja^uj)) is purely real. 

Fig. H shows the FT Z-map Z(q, uj) for a representa- 
tive set of octet q-vectors at three levels of hole doping, 
for an impurity that modulates the superconducting gap 
and bond parameters (panels a-c), or site energy (panels 
d-f). The value of \i in Eq. [2l] is varied to control the 
hole doping, and the superconducting gap maximum is 
set to Ao = 35 meV in all cases. Changing the doping 
in the system moves the van Hove singularity closer to 
the superconducting gap edge (more hole doped) or fur- 
ther away (less hole doped), as shown in Fig. 2^a), which 
affects the particle- hole asymmetry. As a result, for all 
octet q-vectors, the FT Z-map intensity becomes larger 
as the system is more strongly hole-doped. 

The FT Z-map intensity at qi (Fig. [4^), q4, and qs 
(not shown, but qualitatively similar) has low intensity 
until near the energy at which the ends of the CCE cross 
the AFZB, when the intensity rapidly increases. Recall 
that the QPI intensity is extinguished at the AFZB for 
bond modulating impurities, so this increase is solely due 
to the LDOS modulation from superconducting gap vari- 
ation (see Fig. [3]). The FT Z-map intensity at q2, q7 
(Fig. H b-c) , q3, and q6 (not shown, but qualitatively 
similar) , peaks and then decreases with energy due to the 
decrease in a(uj), until the intensity is completely extin- 
guished at the AFZB for symmetry reasons as discussed 
in Section III. 

The behavior of the octet q-vectors in the case of site 
energy modulation (Fig. [4] d-f), is qualitatively similar 
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FIG. 4. FT Z-map intensity for three different levels of hole doping. The QPI intensity in panels (a)-(c) is extracted from 
a Z-map calculation with a "patch" impurity that modulates the bond and superconducting gap parameters. The energy at 
which the ends of the CCE cross the AFZB for each doping level are denoted by dashed lines. The QPI intensity in panels 
(d)-(f) is extracted from a calculation with a "patch" impurity that modulates the site energy The behavior is qualitatively 
similar to that of bond/gap modulation. The QPI intensity is in arbitrary units. 



to that for bond/gap modulation. There is no symmetry 
reason for the extinction of FT Z-map intensity at the 
AFZB, so q2 evolves smoothly through this boundary. 
Because the behavior of the Z-map intensity for site en- 
ergy modulation and bond/gap modulation are similar, 
the rest of the FT Z-map figures shown in the paper are 
computed with only bond/gap modulations. A compar- 
ison of the FT Z-maps for these three types of impurity 
scatterers, all considered separately, has also appeared in 
Ref.EU 

To summarize our results on the energy dependence 
of the FT Z-map intensity, we have shown that the en- 
ergy and doping evolution of the octet q-vectors reflect 
the particle-hole asymmetry in the underlying bandstruc- 
ture. We also have shown that the different types of im- 
purity scatterers give different q-dependent QPI intensi- 
ties, so that the intensity of QPI at certain wavevectors 
could be used to identify the type of underlying impurity 
scattering. Namely, q-vectors that rise quickly from neg- 
ligible intensity at the AFZB can be associated with gap 
modulating impurities, while those that peak and then 
extinguish at the AFZB are from bond modulations; site 
energy modulating impurities can exhibit both of these 
behaviors. While these behaviors offer a way to distin- 
guish between impurity types, disentangling the effects 
of multiple impurity types in one system is still difficult. 



In this case, techniques that highlight certain scattering 
types, such as the use of magnetic field to enhance octet 
q-vectors arising from gap inhomogeneity scattering, as 
demonstrated in Refs.HQland |4T| are especially relevant. 

In order to facilitate the comparison of FT Z-maps 
with DWs in the subsequent sections to the case of pure 
superconductivity discussed here, we plot the FT Z-map 
Z(q, uj) and spectral function A(k, uo) at three bias en- 
ergies in Figure [5] I n all subsequent FT Z-map plots in 
this paper (Figs~|5j [8| [lO|-[T2l), the filling is set to 
< n >=0.73 and the Z(qju) intensity color scale for all 
panels in each figure is the same, although this scale has 
been adjusted between figures for clarity. Note that for 
each point in Fig. |1J the QPI intensity at each octet q- 
vector was extracted from FT Z-maps like these. Peaks 
at wavevectors q2, q3, q6, and q^j are present at low 
energies and then extinguish as energy increases, while 
wavevectors qi, q4 and qs appear only at high energy. 
The peaks all disperse according to the octet model. 

V. LONG RANGE DENSITY WAVE ORDER 

We now explore the signatures of DW order in the 
Z-map. This is motivated by the observation that the 
Z-map intensity is sensitive to particle-hole asymmetry, 
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FIG. 5. d-wave superconductivity. Z(q, cj) (left column) 
and spectral function A(k, lS) (right column) for |cj/Ao| = 
0.075 (a and b), |cj/A = 0.5| (c and d), and \uj/A dw = 0.85| 
(e and f). Here the magnitude of the d-wave superconduct- 
ing gap is Ao = 35 meV. Blue circles mark the locations of 
the octet q- vectors on Z(q, cj), and the ends of the CCE on 
A(k, lj), respectively. 



as discussed in the previous sections in the context of a 
particle-hole asymmetric bandstructure. Since DWs also 
break particle-hole symmetry, the Z-map may be a useful 
quantity to consider in DW systems. The Z-map could 
be used to study any material system with a charge or 
spin order; for the case of the cuprates, it could be used 
to search for an underlying DW order in the pseudogap 
phase. 

Here we consider charge order, although spin orders 
could be treated in an analogous manner. We study two 
different DW orders in this work, illustrated in Fig. |6j 
commensurate Qi = (tt, tt) order, and incommensurate 
checkerboard order. The wavevectors of the checkerboard 
order, nesting the fiat regions of the Fermi surface near 
the antinode for the bandstructure defined in Eq. [2TJ are 
Qf = (0.26tt,0), Q| = (0,0.26tt)P We focus on these 
two DW orders because the (tt, tt) case is the simplest to 
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FIG. 6. DW ordering vectors: commensurate Qi = (tt, tt) 
(red) and checkerboard Q| = (0.26tt,0), Q| = (0, 0.26tt) 
(blue). 



treat, making the underlying physics most clear, while 
the checkerboard order has been proposed as a possi- 
ble explanation for the checkerboard modulations in the 
cuprates. The goal of the next few sections is to identify 
the signatures of these DW orders in the Z-map, and de- 
termine how to differentiate superconducting QPI from 
them. In this section we focus on long range DW order, 
while in subsequent sections we consider the coexistence 
of superconductivity with DWs, as well as fluctuating 
DW order. 

We first consider a Q = (tt, tt) DW with an isotropic 
gap Ajjw The Hamiltonian is expressed in the Nambu 
formalism (Eq. |8|, where now ip^o- = (ck<r? C k+Qcr) and 

A / e k A DW \ 

and the sum in Eq.[8]runs over the reduced Brillouin zone 
\k x \ + \k y \ < tt. While the mean- field treatment of the 
DW order used in this paper certainly cannot fully ex- 
plain the pseudogap in the underdoped cuprates, where 
electronic correlations, neglected in this work, are impor- 
tant, it can offer insight into generic signatures of DW 
order. The CCE are ellipses, as depicted in Fig. [7J Be- 
cause a DW order breaks particle-hole symmetry, unlike 
superconductivity, the CCE at bias energies +u and — u 
are of different sizes. Since the CCE are not as sharply 
pointed as the superconducting banana-shaped CCE, the 




FIG. 7. CCE for Q = (tt.tt) long range DW order with 
an isotropic gap A^v^=100 meV. The normal state Fermi 
surface is shown with dashed lines. 



maxima in Z(q, uj) are more arc-like, as discussed below. 

We now turn to the signature of a Q = (71", tt) DW 
in the FT Z-map, which is computed with the same 
self-consistent T-matrix formalism described in Section 
IV, except now the Green's function is derived from the 
Hamiltonian given in this section. While in the supercon- 
ducting state, maxima in n(r,cj) correspond to minima 
in n(r, — cj), in the DW state the modulation patterns in 
n(r, uj) and n(r, — uj) are fundamentally different. As a re- 
sult, Z(r,o;) does not enhance the strength of peaks as it 
does in the superconducting case, rather, separate peaks 
appear corresponding to maxima (minima) in n(r, ±gj). 

The FT Z-map for a Q = (tt, tt) DW with gap A DW = 
100 meV is shown in Fig. [8] (panels a, c, and e) for 
three bias energies. The sum of the spectral functions 
at positive and negative bias energy, A(k, uo) + A(k, —a;), 
is displayed in panels b, d, and f. By considering this 
sum, the difference between the spectral function at pos- 
itive and negative bias energy is clear. A "double- arc" 
structure appears in Z(q, gj) at all energies, and becomes 
more pronounced with increasing bias energy. This struc- 
ture arises because the CCE at positive and negative 
bias energy are of different sizes, and the size difference 
becomes more pronounced with increasing energy. The 
color coded arrows in Fig. [8] mark areas of high inten- 
sity in Z(q,cj), and the corresponding scattering vectors 
connecting regions of high intensity in the spectral func- 
tion from which they arise. At low energy (Fig. [8^1, b), 
scattering is dominated by a wavevector that connects 
the ends of the ellipses and a wavevector that connects 
the ellipse centers. Note that because the CCE at ±cj 
are slightly different, there is a smearing of intensity in 
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FIG. 8. Long range Q = (tt, tt) DW. Z(q, uj) (left column) 
and sum of spectral function at positive negative bias energies, 
A(k, uj) + A(k, — uj) (right column) for \uj/Adw\ =0.1 (a and 
b), \uj/Adw\ = 0.3 (c and d), and \uj/Adw\ = 0.6 (e and 
f). Here Adw = 100 meV. Peaks in Z(q, cj) are marked with 
arrows color coded with the corresponding scattering vector 
in the spectral function. 



Z(q,cj) (at the blue arrow). The same two wavevectors 
dominate as energy increases; however, at intermediate 
energy (Fig. ^p, d), the single peak marked by the blue 
arrow in Fig7^8^i now splits in two (blue and purple ar- 
rows), because the nesting vectors for the CCE at ±u are 
now significantly different, as shown in Fig. [HH. This ef- 
fect is enhanced with increasing energy (FigT|8e, f). The 
appearance of these "doubled" peaks in the Z-map is an 
indication of DW order, rather than superconductivity. 
We note that a doubling of the number of QPI octet 
wavevectors from particle-hole asymmetry has also been 
discussed in Ref. |25j 

For the case of checkerboard order, the Fermi surface 
is significantly reconstructed due to the many higher or- 
der harmonics associated with an incommensurate order. 
Again, due to particle-hole symmetry breaking in the the 
DW phase, the CCE at ±co> are different, so the "dou- 
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bling" of peaks in the FT Z-map, discussed for the (ir, tt) 
DW, still holds. However, many peaks would appear in 
the FT Z-map due to the significant Fermi surface re- 
construction, so it may be difficult to differentiate the 
"doubled" peaks. 



VI. COEXISTENCE OF 
SUPERCONDUCTIVITY AND LONG RANGE 
DENSITY WAVE ORDER 

We now consider superconductivity coexisting with 
long range Q = (7r,7r) DW order. Since the Z-map en- 
hances the superconducting QPI peaks due to the co- 
herence factors, while it creates a "doubling" of features 
due to DW order, a system with both phases will show 
an interplay of these two effects. Again, while here we 
focus on the particular case of the cuprates, the Z-map is 
potentially useful for studying any material system with 
possible DW orders coexisting or competing with super- 
conductivity. 

The Nambu spinors in the Hamiltonian defined 
in Eq. [i] now have four components: ^ = 

( c kf c k+Qt' c - k i' c - k -Qi)' The matrix Ak is 4 x 4: 



/ Adw Ak 

A k -e k 
\ -A k A DW 




where the sum in the Hamiltonian is over the reduced 
Brillouin zone 1^1 + 1^1 < tt. The CCE for coexisting su- 
perconductivity and (tt, tt) DW order are shown in Fig. [9] 
Due to the unit cell doubling in the DW phase, shadow 
CCE appear outside the reduced Brillouin zone, although 
these are suppressed in the spectral function due to the 
coherence factors. Because Q = {tt, tt) does not nest 
the nodal region, at low energies, the CCE are banana- 
shaped as in the pure superconducting state. Upon in- 
creasing energy, the banana-shaped CCE evolve into el- 
lipses reminiscent of the CCE in the pure DW phase. 

The FT Z-map is evaluated using the self-consistent 
T-matrix formalism described in Section IV, except the 
T-matrix and Green's function are now 4x4. The FT 
Z-map and the sum of spectral functions at positive and 
negative bias energies, for coexisting superconductivity 



and a Q = (7r, 7r) long range DW are shown in Fig. 10 
The DW has an isotropic gap A^^=100 meV, while the 
(i-wave supeconducting gap maximum is Aq =35 meV. 
At low energy the DW order has little effect on the spec- 
tral function (compare Fig. 10 3 and Fig. [5Jd) , because 
the DW vector does not nest the nodal region. As en- 
ergy increases and the DW order causes a bend back 
in the Fermi surface, the spectral function at positive 



and negative bias energies differs (Fig. |10f). The ends of 
the banana-shaped superconducting CCE from Fig. [5] are 
marked by blue dots, note that the region of maximum 
spectral intensity now shifts away from these points. 




FIG. 9. CCE for coexisting superconductivity and 
Q = (7r, 7r) long range DW with A = 200 meV and 

Adw — 250 meV (values chosen for clarity). The normal 
state Fermi surface is shown with dashed lines. 



Reflecting the energy evolution of the spectral func- 
tion, at low energy the FT Z-map is fairly similar to 
the pure superconducting case (compare Fig. 10 1 and 
Fig. [5^i), although a strong arc-like feature appears due 
to scattering between the centers of the CCE (green ar- 
row). With increasing energy, the FT Z-map deviates 
more from the superconducting case, so that in Fig. [lOfc 
the scattering vectors nesting the ends of the elliptical 
CCE show up (purple arrow). Finally, in Fig. 10 3, the 
FT Z-map is dominated by the scattering vectors due to 
the DW, discussed in the previous section, rather than 
the octet q-vectors. 



VII. COEXISTENCE OF 
SUPERCONDUCTIVITY AND FLUCTUATING 
DENSITY WAVE ORDER 

Despite extensive efforts, the only long range DW or- 
ders in the cuprates identified thus far with bulk sensitive 
probes are antiferromagnetic spin order in the undoped 
insulating compounds, and static spin stripe order in 
doped La( 2 -^- 2 y) (Sr,Ba) x (Nd,Eu) 2y Cu04. 5 ° As a result, if 
the pseudogap arises from a DW order, this order would 
likely be of a fluctuating type, which cannot be picked up 
by these probes. Motivated by this, we now consider the 
scenario of superconductivity coexisting with fluctuating 
DW order. We again consider two types of DW order: 
commensurate Q = (tt, tt) order, and incommensurate 
checkerboard order with wavevectors Qf = (0.267T, 0), 
Q! = (0,0.26tt). 
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FIG. 10. Coexisting d-wave superconductivity and 
Q = (7r,7r) long range DW order. Z(q, a;) (left column) 
and sum of spectral function at positive negative bias energies, 
A(k,u) + A(k, -a;) (right column) for |cj/A | = 0.075 (a and 
b), \oj/A \ = 0.5 (c and d), and \oj/A \ = 0.85 (e and f). Here 
Ao = 35 meV and Adw — 100 meV. Blue circles mark the 
locations of the octet q- vectors on Z(q, cj), and the ends of the 
CCE on A(k, uj), respectively. New features arising from the 
DW order in Z(q, uj) are marked with arrows color coded with 
the corresponding scattering vectors in the spectral function. 



FIG. 11. Coexisting d-wave superconductivity and 
Q = (7r,7r) fluctuating DW order. Z(q, uj) (left column) 
and sum of spectral function at positive negative bias energies, 
A(k,uj) + A(k, -uj) (right column) for |cj/A | = 0.075 (a and 
b), |cj/Ao| = 0.5 (c and d), and |cj/A | = 0.85 (e and f). 
Here Ao = 35 meV and Adw = 100. Blue circles mark the 
locations of the octet q- vectors on Z(q, cj), and the ends of the 
CCE on /A(k,cj), respectively. New features arising from the 
DW order in Z(q, uj) are marked with arrows color coded with 
the corresponding scattering vectors in the spectral function. 



The fluctuating DW is included as a self-energy term 
in the normal state non-interacting Green's function as 
in Refs.EUandE! 



1 



iuj n r k Ek(iu) n ) 
where the self-energy term is given by 



(27) 



Sk(i^n) = J dqP(q)A 2 DW /(iuj n - e k+q ) (28) 

The function P(q) is a Lorentzian that peaks at the DW 
ordering vectors. The superconducting Green's function 



is now given by 



G 1 (k,iu n ) = iu n I- 



e k + ^k(^n) 
A k 



A k 

-e k - H k (-icj n ) 



This modified Green's function is now inserted into the 
self-consistent T-matrix formalism described in Section 
IV. 

First consider the case of a fluctuating (7r,7r) DW co- 
existing with superconductivity (Fig. [TT| ) . Compared to 
long range order (Fig. 10), fluctuating order has less ef- 



fect on the FT Z-map, at least at low energy. The spec- 
tral function is similar to the spectral function without 
DW order (Fig. [XT})), whereas at higher energies there 
is a significant Fermi surface bend back and redistribu- 
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tion of spectral weight (Fig. [TTJl, f). The FT Z-map at 
low energy (Fig. [TT^ i) is nearly identical to the FT Z- 
map with superconductivity only (Fig. [5|i), except the 
intensity at the octet q-vectors is already reduced, re- 
flecting redistribution of spectral weight by the DW. In 
Fig. [TTfc , new features start to appear in the FT Z-map, 
and finally in Fig. [TTfc the "doubled" peaks from the DW 
order completely dominate the FT Z-map (purple and 
blue arrows). The signatures of fluctuating (7r,7r) DW 
order coexisting with superconductivity are qualitatively 
similar to those observed for long range order, with su- 
perconductivity dominating at low energies and the DW 
order dominating at larger energies. 

The fluctuating checkerboard order coexisting with su- 
perconductivity (Fig. 12) causes a significant breakup of 
the spectral function at all energies, which is plotted sep- 
arately for positive and negative bias energy for clarity. 
Note that the spectral function at ±cj looks increasingly 
different as \u\ grows. Intensity at the octet q-vectors is 
visible in Fig.[l2^, however the FT Z-map in Fig. 12 1 and 
g is completely dominated by wavevectors introduced by 
the shifted spectral weight. Due to the differences be- 
tween A(k, uS) and A(k, — cj), the Z(q, cj) intensity has 
separate contributions from each. However, the large 
number of scattering vectors obscures the "doubling" ef- 
fect clearly visible for (tt, tt) ordering. 

To summarize, these results show that, like a long 
range DW order, a fluctuating DW order has a signif- 
icant effect on the FT Z-map, and different DW orders 
produce markedly different results. Due to particle-hole 
symmetry breaking in the DW phase, a "doubling" of 
peaks occurs in the FT Z-map, which is evident for a 
simple (7r, tt) order. The Fermi surface is significantly 
reconstructed in the case of a checkerboard order, intro- 
ducing a large number of new peaks into the FT Z-map. 
As a result, the FT Z-map is a useful tool to both iden- 
tify the presence of a DW order (for the case of simple 
orders, where the peak "doubling" is clearly visible), and 
to distinguish between different types of orders. 



VIII. CONCLUSIONS 

In this paper we presented a systematic study of the 
signatures of superconductivity and DW order in the Z- 
map. For superconducting QPI, we showed how impuri- 
ties that modulate the bond, superconducting gap, and 
site energy parameters can be differentiated by the differ- 
ent QPI patterns they produce. We noted that due to its 
definition, the Z-map is inherently sensitive to particle- 
hole asymmetry, and showed that the evolution of its in- 
tensity with energy reflects the underlying particle-hole 
asymmetry of the cuprate bandstructure. 

In the second part of the paper, motivated by recent 
experimental results suggesting a possible DW origin for 
the pseudogap in the cuprates, we explored the effect of 
both long range and fluctuating DW orders on the FT Z- 
map. The reorganization of the Fermi surface by the DW 



order introduced new peaks into the FT Z-map. Due to 
the particle-hole symmetry breaking in a DW phase, scat- 
tering wavevectors at positive and negative bias energies 
are different, leading to a "doubling" of the wavevectors 
in the FT Z-map. For a simple (7r,7r) order, the features 
can be easily connected to scattering vectors between re- 
gions of large intensity in the spectral function. For the 
more complicated checkerboard order, this is not possi- 
ble. However, we note that different types of DW order 
produce very different signatures, providing a means of 
distinguishing between them. 

As is evident from this discussion, the Z-map contains 
a great deal of information. However, due to this com- 
plexity, features due to different types of impurity scat- 
tering and particle-hole asymmetries of various origins 
may obscure the signals from underlying DW orders. In 
addition, STS resolution issues and tunneling matrix el- 
ements, not considered here, may complicate detection 
of the subtle differences in the Z-map patterns from su- 
perconductivity and different types of DW order. This 
could explain why the signatures discussed in this paper 
have not been detected, even if present in the system. 

Finally, we note that in the phase diagram of many 
materials, there is a proximity between DW order and 
superconductivity, for example, the transition from anti- 
ferromagnetic spin order to superconductivity as a func- 
tion of doping in both the cuprate and pnictide families 
of superconductors, coexisting CDW order and conven- 
tional superconductivity in NbSe2, 53 and the pressure- 
induced transition from a CDW to a superconducting 
state in the rare earth tritelluride compounds (RTes)}^ 
As a quantity that can distinguish between these phases, 
the Z-map may be useful in unravelling the relation be- 
tween superconductivity and charge and spin orders in 
many compounds. 
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FIG. 12. Coexisting d-wave superconductivity and fluctuating checkerboard DW order. Z(q, uj) (left column), 
spectral function at positive bias energy A(k, uj) (middle column), and spectral function at negative bias energy, A(k, — uj) 
(right column) for |cj/A | = 0.075 (a-c), \uj/A \ = 0.5 (d-f), and |cj/A | = 0.85 (g-i). Here A = 35 meV and Adw — 100 
meV. Blue circles mark the locations of the octet q- vectors on Z(q, cj), and the ends of the CCE on A(k, uj), respectively. New 
features arising from the DW order in Z(q, uj) are marked with arrows color coded with the corresponding scattering vector in 
the spectral function. The checkerboard wavevectors are set to Q| = (0.267r,0) and Q| = (0,0.267r). 



1 0. Fischer, M. Kugler, I. Maggio-Aprile, C. Berthod, and 
C. Renner, Reviews of Modern Physics 79, 353 (2007). 

2 M. Randeria, R. Sensarma, N. Trivedi, and F.-C. Zhang, 
Phys. Rev. Lett. 95, 137001 (2005). 

3 P. W. Anderson and N. P. Ong, Journal of Physics and 
Chemistry of Solids 67, 1 (2006). 

4 M. R. Norman, D. Pines, and C. Kallin, Advances in 
Physics 54, 715 (2005). 

5 A. Damascelli, Z. Hussain, and Z.-X. Shen, Rev. Mod. 
Phys. 75, 473 (2003). 

6 M. Hashimoto, R.-H. He, K. Tanaka, J.-P. Testaud, 
W. Meevasana, R. G. Moore, D. Lu, H. Yao, Y. Yoshida, 
H. Eisaki, et al., Nature Physics 6, 414 (2010). 



7 R.-H. He, M. Hashimoto, H. Karapetyan, J. D. Koralek, 
J. P. Hinton, J. P. Testaud, V. Nathan, Y. Yoshida, H. Yao, 
K. Tanaka, et al., Science 331, 1579 (2011). 

8 A. V. Balatsky, I. Vekhter, and J.-X. Zhu, Rev. Mod. Phys. 
78, 373 (2006). 

9 H. Alloul, J. Bobroff, M. Gabay, and P. J. Hirschfeld, Rev. 
Mod. Phys. 81, 45 (2009). 

10 C. Howald, P. Fournier, and A. Kapitulnik, Phys. Rev. B 
64, 100504 (2001). 

11 S. H. Pan, J. P. O'Neal, R. L. Badzey, C. Chamon, H. Ding, 
J. R. Engelbrecht, Z. Wang, H. Eisaki, S. Uchida, A. K. 
Gupta, et al., Nature (London) 413, 282 (2001). 

12 J. E. Hoffman, K. McElroy, D.-H. Lee, K. M. Lang, 
H. Eisaki, S. Uchida, and J. C. Davis, Science 297, 1148 



14 



(2002). 

13 J. E. Hoffman, E. W. Hudson, K. M. Lang, V. Madhavan, 
H. Eisaki, S. Uchida, and J. C. Davis, Science 295, 466 

(2002) . 

14 K. M. Lang, V. Madhavan, J. E. Hoffman, E. W. Hudson, 
H. Eisaki, S. Uchida, and J. C. Davis, Nature (London) 
415, 412 (2002). 

15 K. McElroy, R. W. Simmonds, J. E. Hoffman, D.-H. Lee, 
J. Orenstein, H. Eisaki, S. Uchida, and J. C. Davis, Nature 
(London) 422, 592 (2003). 

16 Y. Kohsaka, K. Iwaya, S. Satow, T. Hanaguri, M. Azuma, 
M. Takano, and H. Takagi, Physical Review Letters 93, 
097004 (2004). 

17 T. Hanaguri, C. Lupien, Y. Kohsaka, D.-H. Lee, 
M. Azuma, M. Takano, H. Takagi, and J. C. Davis, Nature 
(London) 430, 1001 (2004). 

18 M. Vershinin, S. Misra, S. Ono, Y. Abe, Y. Ando, and 
A. Yazdani, Science 303, 1995 (2004). 

19 K. McElroy, D.-H. Lee, J. E. Hoffman, K. M. Lang, J. Lee, 
E. W. Hudson, H. Eisaki, S. Uchida, and J. C. Davis, Phys- 
ical Review Letters 94, 197005 (2005). 

20 T. Hanaguri, Y. Kohsaka, J. C. Davis, C. Lupien, I. Ya- 
mada, M. Azuma, M. Takano, K. Ohishi, M. Ono, and 
H. Takagi, Nature Physics 3, 865 (2007). 

21 Y. Kohsaka, C. Taylor, K. Fujita, A. Schmidt, C. Lupien, 
T. Hanaguri, M. Azuma, M. Takano, H. Eisaki, H. Takagi, 
et al., Science 315, 1380 (2007). 

22 Y. Kohsaka, C. Taylor, P. Wahl, A. Schmidt, J. Lee, K. Fu- 
jita, J. W. Alldredge, K. McElroy, J. Lee, H. Eisaki, et al., 
Nature (London) 454, 1072 (2008). 

23 W. D. Wise, M. C. Boyer, K. Chatterjee, T. Kondo, 
T. Takeuchi, H. Ikuta, Y. Wang, and E. W. Hudson, Na- 
ture Physics 4, 696 (2008). 

24 W. D. Wise, K. Chatterjee, M. C. Boyer, T. Kondo, 
T. Takeuchi, H. Ikuta, Z. Xu, J. Wen, G. D. Gu, Y. Wang, 
et al., Nature Physics 5, 213 (2009). 

25 J. Lee, K. Fujita, A. R. Schmidt, C. K. Kim, H. Eisaki, 
S. Uchida, and J. C. Davis, Science 325, 1099 (2009). 

26 C. V. Parker, P. Aynajian, E. H. da Silva Neto, A. Pushp, 
S. Ono, J. Wen, Z. Xu, G. Gu, and A. Yazdani, Nature 
(London) 468, 677 (2010). 

27 K. Fujita, A. R. Schmidt, E.-A. Kim, M. J. Lawler, D. H. 
Lee, J. C. Davis, H. Eisaki, and S.-i. Uchida, Journal of 
the Physical Society of Japan 81, 011005 (2012). 

28 K. Fujita, I. Grigorenko, J. Lee, W. Wang, J. X. Zhu, J. C. 
Davis, H. Eisaki, S. Uchida, and A. V. Balatsky, Phys. 
Rev. B 78, 054510 (2008). 

29 Q.-H. Wang and D.-H. Lee, Phys. Rev. B 67, 020511 

(2003) . 

30 L. Capriotti, D. Scalapino, and R. Sedgewick, Phys. Rev. 
B 68, 014508 (2003). 



31 I. M. Vishik, E. A. Nowadnick, W. S. Lee, Z. X. Shen, 
B. Moritz, T. P. Devereaux, K. Tanaka, T. Sasagawa, and 
T. Fujii, Nature Physics 5, 718 (2009). 

32 B. M. Andersen and P. J. Hirschfeld, Phys. Rev. B 79, 
144515 (2009). 

33 S. Chakravarty, R. B. Laughlin, D. K. Morr, and C. Nayak, 
Phys. Rev. B 63, 094503 (2001). 

34 S. A. Kivelson, I. P. Bindloss, E. Fradkin, V. Oganesyan, 
J. M. Tranquada, A. Kapitulnik, and C. Howald, Rev. 
Mod. Phys. 75, 1201 (2003). 

35 E.-A. Kim and M. J. Lawler, ArXiv e-prints (2012), 
1205.6408. 

36 J.-X. Li, C.-Q. Wu, and D.-H. Lee, Phys. Rev. B 74, 
184515 (2006). 

37 C. Li, S. Zhou, and Z. Wang, Phys. Rev. B 73, 060501 
(2006). 

38 G. Seibold, M. Grilli, and J. Lorenzana, Phys. Rev. Lett. 
103, 217005 (2009). 

39 D. Podolsky, E. Demler, K. Damle, and B. I. Halperin, 
Phys. Rev. B 67, 094514 (2003). 

40 T. Pereg-Barnea and M. Franz, Phys. Rev. B 78, 020509 

(2008) . 

41 T. Hanaguri, Y. Kohsaka, M. Ono, M. Maltseva, P. Cole- 
man, I. Yamada, M. Azuma, M. Takano, K. Ohishi, and 
H. Takagi, Science 323, 923 (2009), 0903.1674. 

42 D. Zhang and C. S. Ting, Phys. Rev. B 67, 100506 (2003). 

43 D. Zhang and C. S. Ting, Phys. Rev. B 69, 012501 (2004). 

44 T. Pereg-Barnea and M. Franz, International Journal of 
Modern Physics B 19, 731 (2005). 

45 L. Zhu, W. A. Atkinson, and P. J. Hirschfeld, Phys. Rev. 
B 69, 060503 (2004). 

46 T. S. Nunner, W. Chen, B. M. Andersen, A. Melikyan, and 
P. J. Hirschfeld, Phys. Rev. B 73, 104511 (2006). 

47 L. Dell'Anna, J. Lorenzana, M. Capone, C. Castellani, and 
M. Grilli, Phys. Rev. B 71, 064518 (2005). 

48 H.-D. Chen, O. Vafek, A. Yazdani, and S.-C. Zhang, Phys- 
ical Review Letters 93, 187002 (2004). 

49 A. Fang, L. Capriotti, D. Scalapino, S. Kivelson, 
N. Kaneko, M. Greven, and A. Kapitulnik, Phys. Rev. 
Lett. 96, 017007 (2006). 

50 J. M. Tranquada, B. J. Sternlieb, J. D. Axe, Y. Nakamura, 
and S. Uchida, Nature (London) 375, 561 (1995). 

51 P. A. Lee, T. M. Rice, and P. W. Anderson, Physical Re- 
view Letters 31, 462 (1973). 

52 N. Harrison, R. D. McDonald, and J. Singleton, Phys. Rev. 
Lett. 99, 206406 (2007). 

53 S. V. Borisenko, A. A. Kordyuk, V. B. Zabolotnyy, D. S. 
Inosov, D. Evtushinsky, B. Biichner, A. N. Yaresko, 
A. Varykhalov, R. Follath, W. Eberhardt, et al., Phys. 
Rev. Lett. 102, 166402 (2009). 

54 J. J. Hamlin, D. A. Zocco, T. A. Sayles, M. B. Maple, 
J. H. Chu, and I. R. Fisher, Phys. Rev. Lett. 102, 177002 

(2009) . 



